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0> ■ ABSTRACT 

in ■ 

' We consider the tidal interaction of a binary consisting of a fully convective 

. primary star and a relatively compact mass. Using a normal mode decom- 

o ■ 

■ position we calculate the evolution of the primary angular velocity and orbit 

C^h' for arbitrary eccentricity. The dissipation acting on the tidal perturbation is 

O ' assumed to result from the action of convective turbulence the effects of which 

, are assumed to act through an effective viscosity. A novel feature of the work 

• • . presented here is that, in order to take account of the fact that there is a 

relaxation time, tc, being the turn-over time of convective eddies, associated 



' with the process, this is allowed to act non locally in time producing a de- 

pendence of the dissipation on tidal forcing frequency. Results are expressed 
in terms of the Fourier coefficients of the tidal potential assumed periodic in 
time. We find useful analytical approximations for these valid for sufficiently 
large values of eccentricity e > 0.2. 

We show that in the framework of the equilibrium tide approximation, 
when the dissipative response is frequency independent, our results are equiv- 
alent to those obtained under the often used assumption of a constant time 
lag between tidal response and forcing. 

We go on to consider the case when the frequency dependence of the dis- 
sipative response is cx l/(l + (cij„i_feic)''), where ujm,k is the apparent frequency 
((patternspeed) x m) associated with a particular harmonic of the tidal forcing 
as viewed in the frame corotating with the primary. We concentrate on the 
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case ujm.ktc >> 1 which is thought to bo appropriate to many astrophysical 
apphcations. Wc study numericaUy and analyticaUy the orbital evolution of 
the dynamical system corresponding to different values of the parameter p. 
We present results from which the time to circularize from large eccentricity 
can be found. 

We find that when p < I the evolution is similar to that found under the 
constant time lag assumption. However, the orbital evolution of the system 
with p > 1 differs drastically. In that case the system evolves through a 
sequence of spin-orbit corotation resonances with flr/florb = n/2, where 17^ 
and florb are the rotation and orbital frequencies and n is an integer. When 
p = 2 we find an analytic expression for the evolution of semi-major axis with 
time for arbitrary eccentricity assuming that the moment of inertia of the 
primary is small. 

We confirm the finding of Ivanov & Papaloizou (2004) on the basis of 
an impulsive treatment of orbits of high eccentricity that equilibrium tides 
associated with the fundamental mode of pulsation and dissipative processes 
estimated using the usual mixing length theory of convection seem to be too 
weak to account for the orbital circularization from large eccentricities of 
extra-solar planets. 

Generalisations and limitations of our formalism are discussed. 

Key words: binaries: general - planetary systems: formation - stars: rotation: 
hydrodynamics 

1 INTRODUCTION 

Tidal interaction between orbiting companions produces dissipation and orbital evolution 
towards a circular orbit which is synchronised with the internal rotations. Such a state is 
attained when the system is close enough for the tidal interaction is strong enough (for a 
discussion in the case of close binary stars see Giuricin, Mardirossian & Mezzetti 1984a,b ). 
The recently discovered extra-solar planets with periods less than about 5 days are in near 
circular orbits which is also possibly a result of tidal interaction ( Marcy et al 2001). 

Zahn (1977) isolated two different types of tidal interaction. The first, applicable to 
fully convective stars or planets, described as the equilibrium tidal interaction, applies when 
all relevant normal modes of oscillation of the tidally disturbed object have much higher 
oscillation frequencies than the orbital frequency. In this case dissipation is presumed to 
occur through a turbulent viscosity which would act much like an anomalous Navier Stokes 
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On equilibrium tides in fully convective planets and stars 3 
viscosity. The second interaction described as being through dynamic tides occurs when 
normal modes have oscillation frequencies comparable to the orbital frequency and their 
excitation cannot be neglected. 

The theory of the equilibrium tide with anomalous Navier Stokes viscosity can be related 
to the classical theory of tides based on the idea that the effect of dissipation is to cause the 
tidal response to lag behind the forcing with a constant time lag ( Darwin 1879; Alexander 
1973; Hut 1981). 

The tidal theory based on a constant time lag is often applied to problems where tidal 
interactions are important either in the case of stellar binaries (eg. Lai 1999; Hurley, Tout & 
Pols 2002 and references therein) or of giant planets orbiting close to their central stars (eg. 
Mardhng & Lin 2002). It may be argued that this approach may be apphcable when the 
dissipative process is one that acts as a standard Navier Stokes viscosity ( sec eg. Eggleton, 
Kiseleva & Hut 1998; Ivanov & Papaloizou 2004, hereafter IP, and below). 

However, it has been noted on theoretical grounds that, if the assumed viscosity results 
from turbulence, there is a natural relaxation time associated with it namely the eddy turn 
over time tc- This may be longer than the orbital period in which case the efficiency of the 
dissipation process should be reduced and that there should be a dependence of it on the 
tidal forcing frequency (see eg. Goldreich & Nicholson, 1977; Goldreich & Keeley 1977; Zahn 
1989). Tidal evolution may differ from that found under the constant time lag assumption 
as a result. The importance of this effect for understanding the observed properties of a 
particular binary system has been stressed by Goldman & Mazeh (1994). 

In this paper we investigate this problem in the simplest possible context. We study tidal 
evolution under the assumption that only the fundamental quadrupole modes are important 
for determining the tidal response and that their eigen-frequency is much larger than the 
orbital frequency so that an equilibrium tide approximation can be made. However, we allow 
the dissipation to have a characteristic time scale associated with it and through this produce 
a frequency dependent response in the manner that has been suggested ( eg. Goldreich & 
Keeley 1977) might be appropriate for convective turbulence. We find then that only when 
the dissipative process acts instantaneously are results equivalent to the constant time lag 
assumption obtained. In other situations, we find that there is the possibility of maintained 
spin orbit resonances and accordingly that the state of pseudo-synchronisation during the 
circularization process differs. 

The plan of the paper is as follows. In section 2 we formulate the basic equations for 
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calculating the response to a tidal forcing potential in terms of a normal mode expansion. 
We also formulate the treatment of dissipation using viscosity that incorporates a relaxation 
time, tc, identified with the turn-over time of convective eddies. This allows the dissipative 
process to act non locally in time and become weaker when the tidal forcing frequency as seen 
in a frame corotating with the tidally perturbed object becomes significantly larger than t~^. 
We derive expressions for the rates of change of energy and angular momentum applicable 
to an orbit with arbitrary eccentricity. We introduce the equilibrium tide approximation 
applicable when the only important pulsation mode for the tidal response is the fundamental 
quadrupole mode and the orbital frequency is very much less than the frequency of this mode. 
We show that the results are equivalent to those under the approximation of a constant time 
lag between tidal forcing and response of the perturbed object only in the special case tc — 
which corresponds to dissipation produced through the action of a standard Navier Stokes 
viscosity. 

In section 3 we develop analytic approximations for the Fourier coefficients of the forc- 
ing potential. We consider the orbital evolution calculating the circularization from large 
eccentricities in section 4. We consider the case when the angular momentum associated 
with rotation of the tidally perturbed object is much less than that of the orbit so that 
pseudo-synchronisation is achieved rapidly. We suppose that the frequency dependence 
of the dissipative response is oc 1/(1 + (c<Jm,fcic)^), where uJm,k is the apparent frequency 
((pattcrnspced) x m) associated with a particular harmonic of the tidal forcing viewed in 
a corotating frame, p is a parameter and uJm,ktc » 1- We discuss when is is possible to 
attain spin orbit resonance and the form pseudo-synchronisation takes as p is varied. Both 
numerical and analytic approaches are used. Finally in section 5 we summarise our results 
and discuss them in the context of extra-solar planets. 

2 BASIC EQUATIONS 

For the case of a fully convective planet or star considered in this paper, the theory of tidal 
perturbation is relatively simple. When the body is not rotating, the pressure (p) modes 
do not contribute significantly to the interaction due to the large values of their oscillation 
frequencies and the gravity {g) modes are absent. Therefore, only the fundamental (/) modes 
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On equilibrium tides in fully convective planets and stars 5 

contribute significantly. ^. For a rotating object there arc also incrtial or r modes which 
may give rise to a complicated spectrum and associated eigen-f unctions. There could be 
significant excitation depending on overlap of these with the forcing potential. For simplicity 
this spectrum is neglected in this paper but in principle could be significant when considering 
the theory of the equilibrium tide ( see eg. Papaloizou & Pringle 1978, 1981; Ogilvie & Lin 
2003 for a discussion). 

We also assume that the rotation rate of the star flj. is uniform with axis of rotation 
directed perpendicular to the orbital plane, and is sufficiently small in magnitude when 
compared to a characteristic value of the eigen-frequency cuq associated with the lowest 
order / mode with I — 2,. It can be shown (see, eg. IP) that for our purposes it suffices to 
take into account only the leading and next order terms in the expansion of all variables in 
the small parameter Qr/^o- Also, only the contribution of the leading quadrupole term in 
the expansion of the tidal potential in spherical harmonics and hence modes with / = 2, are 
considered. From now on we call the planet or star on which the tides are raised the primary 
and the perturbing companion the secondary. 

We adopt spherical polar coordinates (r, 9, (f)) with origin at the centre of mass of the 
primary. Under the above assumptions, the Lagrangian displacement of a fluid element of a 
fully convective primary can be written in the form: 

^= E (1) 

m=0,±2 

where ^q" are the solutions of the standard cigen-mode equation for a non-rotating star (eg. 
Tassoul 1978, Christensen-Dalsgaard 1998) 

-u^Uo + C(C) = 0, (2) 

where uq is the eigen-frequency of the fundamental mode in the non-rotating limit and C is 
the standard self-adjoint operator accounting for action of pressure and self-gravity forces 
on perturbations (eg. Chandrasekhar 1964, Linden-Bell & Ostriker 1967). 

It is important to note that since the displacement vector ^ is a real quantity and the 
eigen-functions have a trivial dependence on the azimuthal mode number m and oc e*"**^, 
it follows from equation (1) that the mode amplitudes bm satisfy 

bm = (b-m)*, (3) 

* Note, that there could be also some non-standard modes like eg. the mode appearing in a model of Jupiter with an assumed 
first order phase transition between molecular and metallic Hydrogen (eg. Vorontsov, 1984), etc. 
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where from now on the superscript * denotes the complex conjugate. 
The eigen-functions are orthogonal in the sense of the inner product 

<CIC>= jd^^piiioT-Ql (4) 

where p is the density of the primary star, and are normalised by the standard condition: 

< CIC >= 1- (5) 
They can be presented in the form 

io = i{r)RY2m{e, <t>)er + msirVY^mie, 0)), (6) 

where are spherical harmonics and the components ^{r)^ and ^{r)s are related to each 
other as 

-r^s = • (7) 

ar r 

This condition follows from the fact that hydrodynamical motion induced in a non-rotating 
isentropic star by a forcing potential must be circulation free. 

The equation governing the Lagrangian displacement produced under tidal forcing is (see 
eg. IP) 

^ + 2uo • Vi + C{i) = - VV^' + P, (8) 

where Uq is the rotational velocity of the primary and the viscous force per unit mass is f^. 

The evolution equation for the mode amplitudes hm obtained by projecting onto the 
eigenvectors ^q' is found to be (see IP) 

hm + ^Ihm + 2imPnrbm = + (9) 

Here the dimensionless coefficient (3 determines correction to Uq due to rotation, it has a 
form (eg. Christensen-Dalsgaard 1998 and references therein): 

= r'\^dTp{2iRis + es)- (10) 
Jo 

The forcing amplitude is related to the forcing tidal potential if)' (assumed to be 
oc ex^{im(j))) through 

fl = -j d'xpiCoT ■ W (11) 
For our problem, the quantity can be written in the form 

£ = ^W-e-'C), (12) 
where the overlap integral Q is expressed as (Press & Teukolsky 1977) 
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On equilibrium tides in fully convective planets and stars 7 

g = 2 r"'drpr\CR + ^Cs). (13) 
Jo 

Here the mass of the secondary star is M, D(t) is the position vector of the secondary 
assumed to orbit in the plane 6 = 7r/2 with D{t) = |D(t)| and $(t) being the associated 
azimuthal angle at some arbitrary moment of time t. For \m\ — 2, Wm — ^Stt/IO and 
Wo — —^jnjb. We find it useful below to represent as multiphcation of a constant 
dimensional factor and a time-dependent dimensionless function: 

f^^WrnC<Pm{t), C = , 0„(t) = ^, (14) 

« D{t) 
where a is the semi-major axis and D{t) = D[t)/a. 

The viscous force amplitudes 

The amplitude associated with the viscous force is given by in equation (9). For the 
standard Navier Stokes equations can be directly calculated (see IP and below). How- 
ever, in this paper we consider a more general type of viscous interaction presumed to be 
induced by convective turbulence and for which the expression for must be modified. 
In particular, as has been pointed out by a number of authors (eg. Goldreich & Nicholson 
1977; Goldreich & Keeley 1977), when the characteristic tidal forcing time is much smaller 
than the characteristic turn-over time of convective eddies, tc, the transport of energy and 
momentum by turbulent motions is ineffective and the effective viscosity coefficient should 
be relatively suppressed. In order to describe the mode energy dissipation arising from the 
turbulence, a time interval of extent t^ centred on the local time t should be considered. 
We would like to take this into account in a simple manner and accordingly introduce a 
non-local form of written in the form 

/oo 
dt'i^it - t'){bm{t') + imQrbmit')), (15) 
-oo 

where the fact that bm+iinflrbm is used means that the dissipation operates on the amplitude 
of the velocity field associated with a particular mode as viewed in the frame co-rotating 
with the primary. We specify a particular form of viscosity kernel 7m(^) below (see equation 
69)), and discuss here only its most general properties. 

Since the tidal forcing time and the turn-over time must be estimated in the co-rotating 
frame, and modes with different m have different apparent frequencies in this frame, 'Jmi't) 
should depend both on m and fir- For general linear perturbations, jmi't) is a complex 
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quantity, but the condition jrn{t) = (7-m(^))* must be valid for self-consistency of equation 
(9). 

In the limit tc ^ the viscous interaction is local in time and jmit) — lo^it), where 
7o is the viscous mode damping rate. It has dimensions of inverse time for a normalised 
eigen-f unction. In that case which corresponds to standard Navier Stokes viscosity it may 
be explicitly expressed in terms of the components ^(r)^ and ^{r)s and the density p and 
kinematic viscosity u of the star (IP): 

70 = 4 J^'''" r'drpui^iCn - ^ + 3^)^ + ^((^^ - S^sf + 2^ + m - Csf)}, (16) 
where the prime stands for differentiation with respect to r. 

Energy and Angular Momentum Exchange between primary and orbit 

RecaUing that for the problem on hand only modes with \m\ — 0, 2 are excited it follows 
from equation (9) that when and (|m| = 0, 2) are equal to zero the following quantities 

are conserved: 

£; = m; + | + c.o'(m; + |), (17) 

where we recall bo is real and 

L = 2{i{b2*b2 - 62(^2)*) - 4/3Q,fe2fc2*). (18) 

These quantities represent the sum of the canonical energies and angular momenta associated 
with the different modes (e. g. Friedmann & Schutz, 1978). 

When the effects of dissipation are included, the evolution equations for E and L follow 
from equation (9) as 

E = Et + K, L = Lt + U, (19) 
where 

Ea = b2{m* + b;f^ + bof^, La = 2t{b;f^-b2{f^r), (20) 

where the index a is either T or v. 

The sum of the orbital energy £'0^6, the energy of the oscillation modes E and the energy 
of the primary star Egt is conserved during the evolution of the system as is the sum of 

t For simplicity wc assume hereafter that the tides are raised only in the primary star. The generalisation to the case of two 
tidally interacting stars is straightforward. Also, note that the energy dissipated by tidal friction and finally radiated away 
from the primary is formally considered as a part of Est- 
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On equilibrium tides in fully convective planets and stars 9 
the corresponding angular momenta. If the slow change of orbital parameters due to tidal 
evolution is neglected , the pulsations of the primary star are strictly periodic and the energy 
of the modes E and the associated angular momentum L time averaged over a period doesn't 
grow with time. In this case we have 

{Eorb) = -{Ev), {Lorb) = -{Lu), (21) 

where the angular brackets denote a time average such that for any quantity Q, (Q) — 

^J^or,dtQ{t). 

orb 

Note that the mean rate of energy dissipation {Ediss) differs from {Ei) for a rotating 
primary being given by 

{Ediss)^{K)-nr{L,), (22) 
and must be negative. On the other hand the sign of E^, can be arbitrary. 

Solution for the mode amplitudes 

2. 0. 1 Fourier expansion of the mode and forcing amplitudes 

In order to solve equation (9) for the mode amplitudes we assume that the orbital parameters 
are fixed so that all quantities of interest are periodic functions of time with period Porb- 
Note that a care should be taken with this assumption (see IP and references therein and 
also the discussion), especially for orbits with high eccentricity e ^. 

We use below a Fourier expansion for the mode amplitudes bm and the dimensionless 
tidal forcing amplitudes in the form: 

oo 

E bm,ke'^\ (23) 

fe=— oo 
oo 

0m = E <t>m,ke'^\ (24) 

fc=— oo 

where we introduce the dimensionless time variable r = Vtt, and Vt is the orbital frequency: 



where Mp is the mass of the primary, and q — Mp/M is the mass ratio. 

t For a highly eccentric orbits this assumption can breaJc down when the excitation of so-called dynamic tides is significant (eg. 
Press & Teukolsky 1977). Also the stochastic evolution of orbital parameters (eg. Mardling 1995, IP and references therein) 
may take place. 
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From the condition bm = {b-m)* it follows that 62,^ = (^-2,-fc)*- The amplitude 60 is real 
and therefore 6o,A; — (^o,-fe)*- 

The Fourier coefficients of the dimensionless tidal forcing amplitudes 0^^^ are given by 

1 /■2tt Q—im^{t)-ikT 

<i>m,k = TT TTTTT^ d'T- (26) 



27r Jo D{tf 

We note that D and $ are related to r through D{t) = (1 — ecos(^)) and tan$ = 
Vl — e^sin^/(cos^ — e) with r = ^ — esin(^). 

On account of the fact that D and $ are even function and odd functions of r respectively, 
we note that (j)rn,k is a function only of the eccentricity and is real. We may therefore write 

0m,fe = -^{0im,k — Pm,k), (27) 



where 

1 f^'^ cos(m$(r))cos(A;T) 1 f^'^ sin(m$(r))sm(A;r) , , 

am,k = - / OT -—— , /3m,fe = - / dT . (28 

TT Jo D(r)3 TT Jo D(r)3 

We note obvious properties of the coefficients such as a„i-k = <^m,k, Pm-k = —Pm,k and 
Q!-m,fe = o;^,*;, /3-m,fc = —Pm,k- Furthermore it is important to note that for the case of 
interest with m = 2 and positive values of k we approximately have a2,k ~ P2,k, ^-nd for 
= 0, /3o,fe = 0. 

We find approximate expressions for a2,k and ao,fc for positive values of k, these only 
being required, in the Appendix. We find it convenient below to represent the dimensionless 
tidal forcing amplitudes (p±2 in terms of the real quantities 0_) defined through 

0±2 = 0+Ti0-. (29) 

One can easily check that 0:2,^ {P2,k) are the Fourier cosine (sine) expansion coefficients of 
4>+ {4>-), and Qio,jt are the Fourier cosine expansion coefficients of 0o- 
We introduce the Fourier transform of the kernel function 

dr-fm{r) exp{-ikT) (30) 

-00 

It is convenient to work with a dimensionless form 7^ defined through 

7m,fc = 77m,fc, (31) 

where 7 is some characteristic value of the damping rate. For the case of the standard Navier 
Stokes viscosity, — 0, 7 = 70 and 'ym,k — 1- 

In order to find solutions of equation (9) we first express in terms of the 7^^^ and 
bm,k with the help of equations (15) and (23) so obtaining 
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oo 

Fra = + ^^)lm,khm,ke'^\ (32) 

—oo 

where 

After substitution of equations (23-24) and equation (32) into equation (9), we obtain 

WmC(j)m,k 

bm,k = ^ — -, (34) 

with 

R^u^- {ki^VL + m/?Or - i{k + 2m(7)7^,fc)n. (35) 

We obtain expressions for (£"1^) and {L^,) by substituting (23) and (32) in equations (19 

- 20) and averaging the result over a period 27r in r so obtaining 

00 

{E,) = -2^2 Y^{k{k - 2a)j2,-kb2,-kbl.k + + 2a)j2,kb2,kblk + ^'7o,fe&o,fefeo, J, (36) 
1 

CO 

(L.) = -m Y,{{k - 2a)j2,-kb2,-kb*2,_k - (^ + 2<j)72,fc&2,fc&;, J, (37) 
1 

where we have used the symmetry properties of bm,k and 7m,fc as well as the fact that the 

Fourier coefficients of the potential vanish when A; = and m = 2, to obtain a sum over 

positive integers k only. It is easy to see that 
00 

(Ediss) = E{(^ + '2'^fl2,kb2,kblk + 2af-f2,-kb2,-kbl_k + A;'7o,ik&o,fe&o, J (38) 

1 

is always negative when 7^ ^ > 0. 

Equations (36 - 38) and (34) together with expressions (26) for 0^^^ and the standard 
Keplerian expressions for the orbital energy and angular momentum give a complete set of 
equation for our model. 

The equilibrium tide limit 

2.0.2 General energy and angular momentum exchange rates 

In the approximation scheme corresponding to calculation of the equilibrium (or quasi- 
static) tide, the orbital and rotational frequencies and the mode damping rate are very 
much smaller in magnitude than the oscillation frequency ujq. In this case these can be 
neglected in the denominator R of equation (34). The Fourier mode amplitudes bm,k are 
then simply proportional to 0^^^ such that 

WmC(f)m,k /gnN 
Om,k — 2 ■ V-^yj 
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Substituting (39) in (36) and (37) we obtain 

{K) = -SiE'^, (40) 

(L.) = -S^Ll (41) 
where 

^^=10^' ^^ = T^' ^^^^ 

and 

oo 2 

S, = 2 ^{A;(/c - 2a)%_k<l>l-k + + 2a)72,fc0|, + ^k^lo,k<Plk}^ (43) 
1 

oo 

^2 = 2Y,{{k - 2a)^2,-k<i>l-k -ik + 2a)72,fe0^,J. (44) 
1 

Here we use the definitions of 7 and ^rn,k, the exphcit form of and the fact that 0^^^ 
are real. The quantities and Lq provide characteristic values of the rates of change of 
the energy and angular momentum and set a characteristic time scale of tidal evolution. 
Taking into account the definition of C (equation (14)) and the normalisation condition for 
the eigen-modes (5) one can easily see that these quantities have the correct dimensions of 
energy and angular momentum per unit of time respectively. 

As discussed above, \(p2,k\ ^ I^'^2,-A;| and in a good approximation we can neglect the 
terms proportional to in equations (40-44) and rewrite them in a simpler form: 

00 r) 

{K) ^ -2El Y.{k{k - 2a)72,-fc0^,_, + i:k^%,k'i>lk}. (45) 
1 ^ 

00 

{U) ^ -2Ll Y.{k - 2a)j2,-k<Pl-k- (46) 
1 

We use equations (45) and (46) below for our numerical calculations. 
2.0.3 An alternative form 

There is a form of equations (40-44) which is useful for relating the work here to that of 
Hut (1981) in which he studied tidal interactions of general orbits using a constant response 
time lag approximation together with the assumption of equilibrium tides. We write them 
in the form: 

{K)^-E^A'^i-2a^2), (47) 
(L,) = -L0(^2- 2(7*3), (48) 
where 
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1 °° 2 

*i = 2 ^ ^'{72,-fe(a2,ifc + P2,kf + 72,fe(a2,it - P2,kf + gTo.feao.J, (49) 

2 oo 

*2 = ^ E k{l2,-k{(^2,k + P2,kf - l2,k{(^2,k - P2,kf}, (50) 
^ 1 
-j^ oo 

*3 = E{72,-fc(«2,ik + P2,kf + 72,ik(a2,ik - P2,kf}- (51) 
^ 1 

Here we have expressed 02,±fc in terms of a2,k and /92,jk using equation (27). Obviously, we 
have 5"! = *i - 2(7*2 and 5*2 = *2 - 2(7*3. 

S.O.^ The limit tc^Q 

The hmit tc corresponds to local behaviour and a Navier Stokes viscosity. Assuming 
in this limit that 7m(t) = 7o^(^), we may adopt '^rn,k = 1 in equations (49-51). In this case 
the infinite series (49-51) can be exactly summed by use of arguments exploiting Parseval's 
theorem. 

Then we directly obtain from equations (49-51) that 

oo -| 

*? = E^'(«2,. + /52%+3<.), (52) 



1 

oo 



*° = 2E^«2,fe/92,ik, (53) 



1 

oo 



*^ = E«. + /32%)- (54) 
1 

Using the definitions of am,k and j3m^k (equation (28)) and also of 0+, 0_ and 0o (equation 
(29), one can easily show that the series (52-54) can be expressed in terms of integrals: 

^l = - r dT{ct>l + ct>'-). (57) 

TT JO 

The evaluation of the integrals (55-56) is straightforward with the result: 

TO 8 31 2 255 4 185 g 25 

*S = ^(l + ye^ + ^e^ + ^e^ + ^e^ (58) 

^^ = ^(l + fe^ + fe^ + ^e«), (59) 
*^ = |(l + 3e^ + ^e^), (60) 



where e = \J{1 — e^). 
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When the expressions (58-60) are inserted into equations (47 - 48) the resulting exchange 
rates of energy and angular momentum are identical to the corresponding expressions ob- 
tained by Hut (1981) provided we make the identification 



where Rp is the primary radius, kn is the apsidal motion constant defined as in Hut (1981), 
and Tff is Hut's typical time of tidal evolution. 

Thus, in the limit of constant (frequency independent) viscosity, our normal mode ap- 
proach to the problem is equivalent to the standard constant time lag approach (Alexander, 
1973, Hut, 1981). However, it has the advantage of enabling us to express the unspecified 
evolution time Th in terms of micro-physical quantities such as the kinematic viscosity u 
and the primary density p and quantities related to the normal mode (see also Eggleton, 
Kiseleva & Hut 1998). Note that equation (61) has been obtained by IP in the limit of highly 
eccentric orbits § (see also equations ( 23 - 25 ) of Kumar, Ao & Quataert 1995). 

3 THE FOURIER EXPANSION OF THE TIDAL FORCE AMPLITUDE 

Here we briefly compare the approximate values for the Fourier coefficients a2,k and ao,k ob- 
tained with our analytic approximations derived in the Appendix with the exact coefficients 
obtained numerically. The various coefficients are plotted in figures 1,2 and 3. In figure 1 we 
compare values of the Fourier coefficients obtained analytically with those obtained numer- 
ically for small e = 0.2 and e = 0.3. One sees that ao.jfc decreases monotonically with k and 
that the analytic and numerically obtained values are in a good agreement for practically 
all values of k. 

The difference is largest and of order of 10 percent at A; = 1 and rapidly decreases 
thereafter with k. The coefficient 0:2, fe has a maximum at kmax = 2 for e = 0.2 and at 
kmax = 3 for e = 0.3. For values of /c > k^ax the analytic and numerically obtained values 
differ by at most 10 percent and approach each other rapidly with increasing k. For values 
of A; < kmax the agreement between the analytic and numerically values is less good. 

In figure 2 the same comparison is given for intermediate values of eccentricity e = 0.4, 
0.5 and 0.6. In this case the agreement between numerical and analj^ic curves is much 

§ Note a misprint in IP. Their coefficients /3i and /32 , and accordingly the numerical factor entering their equation (46) analogous 
to equation (61) must be divided on two. 



Th 



AwjMpiGQf 
5 cu^R^^ 



(61) 



© 1994 RAS, MNRAS 000, 1-37 



On equilibrium tides in fully convective planets and stars 15 



T 1 1 1 1 1 < 1 < 1 1 1 < 1 < 1 1 1 1 1 r 



1 -i, 




Figure 1. The coefficients ao,fc (dashed lines) and 02, fc (solid lines) obtained with the analytic approximations together with 
the same quantities obtained numerically (dotted and dot-dashed lines respectively). Curves of the same type with smaller 
values of a^ f, correspond to e = 0.2 and those with the larger value to e = 0.3. 




Figure 2. As in 1 but for e = 0.4, 0.5 and 0.6. Curves of the same type with increasing values of ajn,k correspond to increasing 
values of e. 

better. The disagreement clearly decreases with increasing e and for e = 0.6 the curves for 
Q;o,fe differ by less than or of the order of ~ 1 percent. In the case of a2,k the agreement is 
less good, however, in contrast with the case of smaller e, the disagreement remains small 
even for k < kmax provided that k is sufficiently large. For example the curves corresponding 
to e = 0.6 differ by less than ~ 10 percent for k > 4. 

In figure 3 we show the comparison for e = 0.7, 0.8 and 0.9. The disagreement between 
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Figure 3. As in 1 and 2 but for e = 0.7, 0.8 and 0.9. Curves of tlie same type with increasing values of a^^fc correspond to 
increasing values of e. 

analytic and numerically obtained values is very small for all interesting values of k. The 
numerical and analytic curves corresponding to Q;o,fc practically coincide with each other. 
The numerical and analytic curves corresponding to a2,k and e = 0.9 differ by less than ~ 1 
percent for k > 10. 

In general the analytic and numerical curves are in better agreement for ao^k than for 
a2,k- This occurs because in order to obtain our analytic estimate for a2,k we had to expand 
cos(2$) in powers of the eccentric anomaly ^. We expanded cos(2$) to fourth order in ^. It 
seems that one can obtain a better approximation by expanding to higher order in of ^. 

We point out that our approximate analytic expressions for the Fourier coefficients as- 
sociated with the tidal amplitudes may be useful in a much more general situation than 
the one considered here but where an analytic form of the tidal potential is required and a 
perturber moves on an eccentric orbit. 

Such situations could naturally arise, say, in the problem of the interaction of an accretion 
disc with a perturbing planet or star, or in a problem of the interaction of a super-massive 
black hole binary with a stellar cluster. 
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4 ORBITAL EVOLUTION 

The evolution equations for the orbital parameters follow from equations (21) and equations 
(47), (48) giving two equations for the rate of change of orbital semi-major axis and specific 
angular momentum in the form 

2E0a%T o TN 37r(GM)2g2^(l + g) 1 



GMm 5 MpUjQ a' 



where Lo^b is the specific orbital angular momentum and we recall that a is the semi-major 
axis, Mp is the mass of the primary, q — Mp/M. 

To determine the ratio of the primary rotation rate and the orbital mean motion, a — 
Qr/Q, and hence close the set of equations (62) and (63) we use the law of conservation of 
angular momentum in the form: 

Lorb + I^r — Lin, (64) 

where / is the moment of inertia of the primary per unit mass and Lj„ taken to be the 'initial' 
value of the specific orbital angular momentum at the beginning of the orbital evolution due 
to equilibrium tides. Note that in addition to a, the quantities ^'i, ^'2 and depend on 
eccentricity e which itself may be expressed in terms of a and Lg^b by use of the standard 
relation e — \l — 



GMo" 

It is convenient to rewrite equations (62) and (63) in dimensionless form by use of a length 
scale , ao, which is defined to be the semi- major axis corresponding to the tidal equilibrium 

state of the dynamical system found from the conditions (da/dt = dL^j-b/dt = 0). 

When a = Oo is determined in this way wc also have e = 0, a = 1 in equilibrium. We 
define the dimensionless semi- major axis, x, to he x — a/ao and the dimensionless angular 



momentum to be y = Lgrb/LQ, where Lq — sf^^^- 

In terms of these variables, equations (62) and (63) take the form 
dr 1 

^ = -pr(*, - 2^*2), (65) 

|-JI^(*-2.*3). (66) 

where we introduce the dimensionless time i —t/t^, with 

^ 5 Mpuj^o 

" 3n{GMQ)^^{l + q) ^ ^ 
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The time t,^ defines a characteristic time scale for tidal evolution of the orbital elements, in 
this case due to equilibrium tides. 



We note that in terms of the dimensionless variables, the eccentricity e = \J{^ ~ u'^/x) 
and the law of conservation of angular momentum (64) takes the form 

y + = Vin^ (68) 

where yin — Lin/ ^/GMoq. We recall that cr = Q^/il and a dimensionless moment of inertia 
/ = 7(1 + q)lal 

As for most cases of astrophysical interest, the angular momentum content of the orbit 
dominates that associated with the rotation of the primary, for simplicity, we consider below 
only the case when the primary effectively has low- inertia assuming that / ^ 1. 

As we see below, in this case, the orbital evolution of the binary is essentially independent 
of the value of /. When the viscosity Fourier transform coefficients 7^^^ are specified, the 
set of equations (65), (66) and (68) becomes complete. 

We note that neither the length scale ao or q enter into the above dimensionless set 
of equations or the associated dimensionless initial conditions making the problem scale 
free. That means that solutions of the tidal problem can be scaled in a self-similar manner 
through adjusting this variable. Note that in order to make such a similarity transformation 
in general, the parameters entering into the specification of the viscosity coefficient may also 
need to be scaled appropriately. 

Explicit form of the viscosity coefficients 

Before solving the equations of tidal evolution we should specify the viscosity Fourier trans- 
form coefficients 7m,fe- As we have mentioned above, it was suggested by Goldreich & Nichol- 
son (1977) and Goldreich & Keeley (1977) that the action of an effective viscosity arising 
from convective turbulence should be suppressed when the typical time scale associated with 
convective eddies is much larger than the time scale associated with the global disturbance 
on which it is presumed to act. To take this effect into account in the simplest manner 
possible, we adopt expressions for 7^ ^ of the form 

^ {i + {Z,ktcyy ^^^^ 

where 

^m,k — + rn^lr\ (70) 



© 1994 RAS, MNRAS 000, 1-37 



On equilibrium tides in fully convective planets and stars 19 
is the apparent frequency associated with the global disturbance as viewed in the frame co- 
rotating with the primary, it being assumed that that is well defined. ^ Here p is a constant 
parameter. 

Goldreich & Nicholson (1977) and Goldreich & Keeley (1977) adopted p = 2. However, 
Zahn(1989) has argued that p ~ 1. Here we shall allow for general values of p. Equation (69) 
can be brought into the standard form (31) with 
~ ^ Tl -^70 



where the dimensionless parameter = tctl^^ determines the magnitude of the suppres 



^0 



sion of the viscosity. We assume below that 1. In this limit, many details of the orbital 

evolution are essentially independent of the particular value of this parameter. 

Effective spin orbit resonances 

A feature arising in the limit 3> 1 is that the action of viscosity is most effective for a 
particular forcing term when the rotational and orbital frequencies are such that {k + ma) — 
0. When this satisfied the potential pattern associated with the Fourier component (m, k) of 
the forcing potential co-rotates with the primary. The relative forcing frequency is thus zero 
and there is no suppression of the viscosity. When ^ 1 viscosity acts strongly only very 
close to the corotation resonances with (k + ma) = 0, occurring when [k = l,2,3....,m = 
—2). Hence, the terminology associated with a resonance is used. 

As we shall see below, the qualitative nature of the orbital evolution depends on whether 
the parameter p exceeds or is less than unity. This is related to the resonant behaviour 
considered above. 

Let us suppose that T* » 1 and we are close to a resonance such that the value of 
the dimensionless rotational frequency a is close to some integer, n, divided by two, thus 
cr ~ n/2. Then it follows from equation (43) and (44) that the quantities Si = (^i — 2(7^2) 
and S'2 = (^2 — 2(7^3) are mainly determined by the terms in the respective summations 
for these expressions with k — n and m — —2 which are resonant in the above sense. Later 
we call these terms the 'resonant terms'. 

We have Si^Sl + S'^'' and 82^81 + , where 

^ We stress that in the case of forced oscillations, the responding global disturbance with particular m and k oscillates with 
the frequency of the corresponding Fourier term in the expansion of the tidal forcing amplitude. 
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^ \ / I zj, — It 

{l + n\n-2a\Px-^p/^)' 




(l + Ti'|n-2a|fa;-3f/2)' 



2TP{n - 2a)02 , 



(72) 



and 5'"'', stand for all other (non- resonant) terms in the series. 

When we have the exact equality, a = n/2, the resonant terms are equal to zero. However, 
their absolute values increase very sharply with an increase of the difference \a — n/2\. For 
example, when ^ 1 and \a — n/2\ ~ <^ 1, the resonant terms are proportional to 
T^~^^ in absolute magnitude. Therefore, for the case p > 1 they can, in principal, be much 
larger than the non-resonant terms. On the other hand, in the case of p < 1 the contribution 
of the resonant terms to Si and S2 is insignificant. We shall see below that in the case p > 1, 
the orbital evolution can proceed through a stage where the condition cr ~ n/2 is maintained 
for a long time. Thus, the system can evolve in a state that maintains a specific spin-orbit 
corotation resonance. 

The case p < 1 

When the parameter p is smaller than 1, the orbital evolution of a low inertia primary star 
is, in many details, similar to the case of assumed constant time lag between response and 
forcing potential considered by Alexander (1973) and Hut (1981). This also corresponds to 
the standard Navier Stokes viscosity {tc — 0)- If the initial value of a is sufficiently small, 
because of the small inertia of the primary, the system quickly relaxes to a state where the 
orbital angular momentum is approximately conserved and we can assume that y ^ fti 1. 
In this case the eccentricity and the semi-major axis are related to each other as 



Also because of the low inertia of the primary, its rotation rapidly approaches the value 
given by fir ~ Cpsfl, where aps corresponds to a state of so-called 'pseudo-synchronisation' 
(Hut 1981) and is obtained from equation (66) after setting dy/dt — 0. 
In a state of pseudo-synchronisation, apg is given by 



However, note that both ^2 and ^3 are, in general, functions of aps, and therefore equation 
(74) should be considered as a non-linear algebraic equation for apg. It is convenient to 
express its solution in the form aps = aian where an represents the solution of (74) for the 
case of constant time lag or the standard Navier Stokes viscosity {tc = 0). This is given by 




(73) 



(Hut,1981) 
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Figure 4. The quantity ai = CTps/cTH as a function of e. The line (Ti = 1 corresponds to p = 0, the solid, dashed, dotted 
and dot-dashed curves correspond to p = 1, 0.75, 0.5, 0.25, respectively. Note that the curve with p = 1 exhibits an irregular 
behaviour at large eccentricities e > 0.85. This is most likely to be a numerical artifact related to the technical problem of 
finding roots of equation (74) in the case of a very rapidly varying function. 



2^3 



1 (1 + fe^ + f e'l + ^e6) 



(75) 



;3 (1 + 3e2 + |e4) 

e = — e^) and we have made use of equations (59) and (60). 

We show the dependence of ai on e obtained numerically for different values of p in 
figure 4. We see that the difference between apg and an is rather small for sufficiently 
large eccentricities. This comes about, in simple terms, because the inverse stellar rotation 
frequency Q~} = 1/ (o-psQ) at pseudo-synchronisation is always of the order of a characteristic 
time of periastron passage. 

To obtain an approximate equation for the evolution of the semi-major axis we substitute 
equation (74) into equation (65) obtaining 



dx 1 , ^ 



)■ 



(76) 



dt ' *3 

As we have mentioned above, resonant effects are not important when p < 1. That means 
that to a good approximation we may neglect unity in the denominator of the expressions 
(71) in the limit oo. Then the functions ^i, ^2 and ^3 in equation (76) do not depend 
on the parameter T*. Taking into account equation (73) we see that for a given p the right 
hand side of equation (76) is only a function of the semi-major axis x. In that case, the vari- 
ables in equation (76) can be separated and equation (76) can easily integrated by numerical 
means. We show the result of integration of equation (76) in figures 5 and 6 together with 
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Time 



Figure 5. The dependence of a; = a{t)/ao on dimensionless time t. The soUd curves are the solutions of equation (76) and the 
dotted curves represent tlic solutions of the full dynamical system. Different curves of the same type correspond to p = 0, 0.25, 
0.5, 0.75 and 1. Curves having larger values of x at the same time i have larger values of p. The initial value of a is set to zero 

the results of numerical integration of the full set of equations (65), (66) and (68). One can 
see that the approximate semi-analytic approach gives a very good approximation to the 
solutions of the full set of equations. 

Note that for our numerical integration we set e = 0.9 initially and terminated the 
integration when e = 0.2. We use a very large value of T^, = 10^ and a small value of the 
dimensionless moment of inertia / = 10~^. We have checked that other reasonable values 
of these parameters dot not lead to a different evolution of our dynamical system. We have 
used our approximate expressions for the Fourier coefficients of the tidal forcing amplitude in 
numerical integration of equations (65), (66) and (68) and good agreement between numerical 
and semi-analytic calculations implies that use of the approximate analytic expressions for 
the Fourier coefficients is justified. 

In figure 5 we show numerical solutions of our equations setting initial value of the 
rotational parameter equal to zero, and in figure 6 we consider the case of a high am — 
100. 

From these and similar plots the dimensionless time required to approach circularity 
starting from large eccentricity can be estimated. This is typically ~ 0.2. Both cases show 
quite similar late time evolution. However, in the case of = 100 we also get a short initial 
period of time when the value of semi-major axis increases with time. During this period of 
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Figure 6. Same as figure 5, but the initial value of <t is set to 100. 

time the rotational energy of the primary is transferee! to the orbital energy. We discuss this 
effect in more detail later on. 

4.1 The case p > 1 

In this case the evolution is qualitatively different from when p < 1. This is because of the 
effects of spin orbit resonances become important. To see this let us consider the case p — 2 
where a simple analytic approach is possible in detail. 

4.1.1 pseudo-synchronisation 

Assume that the moment of inertia is very small and the rotation of the primary is close 
to a spin-orbit resonance of the order n so that Act = a — n/2 <^ 1. Then just as for the 
case p < 1 the assumption of low primary inertia means we can consider the orbital angular 
momentum to be a conserved quantity and set dy/dt — 0. 

It then follows from equation (66) that 5*2 = ^2 — 2(7^3 — 0. As discussed above we 
divide S2 into resonant and non-resonant parts writing 82 = 82 + 82^ we must then have 
5*2 = —82''. Using the explicit form of ^'2 (see equation (72)) we obtain 

l + 4T2x-3/2(Aa)2-4T20|_„- ^ ^ 

The solution of this equation gives the deviation of a from the value corresponding to 
exact resonance as 



© 1994 RAS, MNRAS 000, 1-37 



24 P. B. Ivanov, J. C. B. Papaloizou 



( qnr\2 



where we do not consider the second unphysical root. For a physically meaningful solution 
the expression contained within the square root in equation (78) must be positive. This 
implies that the dynamical system can stay in the resonance only if 

and therefore 



< x^/'T,, (79) 



^3/2 

\Aa\< Aamax = (80) 

Condition (79) tells that for a given T* and x, evolution of the dynamical system is possible 
only if the ratio is sufficiently small. This can always be satisfied if T* is sufficiently 
large. 

However, for a sufficiently large n this ratio increases with decrease of eccentricity e. 
Therefore, during the process of tidal circularization, a resonance with some particular order 
n — 2a becomes at some point impossible. The dynamical system may then rapidly relax to 
some other resonance with n' < n where the ratio is smaller and the condition (79) is satisfied. 
Thus, the analytic approach indicates that a low inertia primary evolves through a sequence 
of spin-orbit resonances with decreasing n during the process of tidal circularization. As we 
see below this character of the evolution is confirmed by numerical simulation. Note that 
the qualitative character of the tidal evolution remains the same for any p > 1. 

4.. 1.2 Evolution of the semi-major axis 

Now let us discuss the evolution of the semi-major axis. We assume that condition (79) is 
satisfied and the primary evolves in a state of spin-orbit resonance with some particular 
n — 2(7. In order to find an evolution law for the semi-major axis we have to evaluate 
Si — '^i — 2(T^2- In general, this quantity can be represented only in terms of complicated 
summations (see equation 43). However, for the case p — 2 these series can be performed 
in the limit T^, — > oo. The resulting expression is remarkably simple. Then we can go on to 
obtain an analytic expression for the dependence x{t) (see below). 

The quantity 5*1 can also be divided on the resonant and non-resonant parts (see equation 
72): Si — 31 + S'^'^ . It is easy to see that — and from the condition of constant orbital 
angular momentum it follows that S2 — —SY, and we have 
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= S'l' - nSY- (81) 

Now the explicit form of 5"! can be written down with help of equations (43) and (44) in the 
form 

S^ = 2Y,{k-n){k- 2a)72,_fe0|_fe + {k + n){k + 2a)%k4,k + ^^'7o,fe0o,if (82) 

In equation (82) we can set 2a = n. We can also neglect unity in the denominators of the 
expressions for 7^^^ and set 2a = n there (see equations (71)) so as to obtain 

"^''^'^Jk^^' 7o,fc=p, (83) 
Using equation (83) in equation (82) gives 

oo Q oo 1 

S, = 2x'Y.i<t>l-k + <Plk + i^<l>l,k) = + f^lk + :^<k), (84) 

k=i k=i 

where in the last equality we use expressions of (f)m,k in terms of a„t,A: and f3rn,k (see equation 
(27)). Note that the series (84) does not contain any divergence aX k = n and we include 
the term with k — n m. the series, assuming with justification when n is large, that it gives 
a negligible contribution. 

According to the Parseval's theorem we have 

oo 1 „27r 1 1-271 

EKk + f^lk)-^- dT{<l>l + cl>l), 2<o + E<. = -/ dT<l>l. (85) 

k=l ^ ''^ k=l ^ ''^ 

Here we recall that the coefficients 0:2,0 and /32,o vanish and the inclusion of these terms in 
the summations is redundant. Furthermore the effect of ao,o in the second summation is 
negligible for orbits of significant eccentricity. From definitions of the quantities 0± and 0o 
(see equation (29)) we have 0^ + 0^ = 0q. Accordingly, both integrals in equation (85) may 
be taken to be *g = ^(1 + 36^ + fe'') (see equation (60)) with e = ^'(l - e^), and 

Taking into account that in our case the orbital angular momentum is approximately con- 
served we can express e in terms of x using equation (73) and substituting equation (86) in 
equation (65) we obtain a remarkably simple equation for evolution of the semi-major axis 
which takes the form 

|._|.-3. I). (S7) 

Note that equation (87) can be easily generalised to take into account the term q;o,o and 
also the terms determined by a non-zero rotational angular momentum of the primary. The 
solution of equation (87) can be expressed in terms of elementary functions 
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Figure 7. The dependence of a; = a/ao on time t. Different curves correspond to different values of 'initial' rotational 

parameter cri„ set in the beginning of computations. The solid, short-dashed, long-dashed, dot-dashed, and dot-dot-dashed 
curves correspond to cri„ = 0, 10, 50, 75 and 100, respectively. The dotted curves represent the analytic solution (88) with 
different values of the integration constant C. Note that the curves corresponding to iTi„ = and 10 (the solid and short-dashed 
curves, respectively), practically coincide. 

where C* is a constant of integration and 







\/x + ^/x_ 


C — 






\/x- ■s/x_ 





(89) 



where x± — |(1±2^ ^) and c± — x±^/2. It is very important to note that solution (88) does 
not depend on the order n of some particular resonance. Therefore it describes approximately 
not only the evolution of the dynamical system in a state of spin-orbit resonance with 
particular value of n, but also a system evolving from one resonance to another but spending 
most of its time in a state of spin-orbit resonances with different values of n. 



Numerical calculations 

We solve equations (65) and (66) numerically for different initial values of the rotational 
parameter and p = 2, and also for (Tj„ = and p = 1.5 and p = 2.5. The dimensionless 
parameter is set to 500 and the dimensionless moment of inertia / is set to 1/300 for 
all numerical calculations. We have checked that other values of these parameters do not 
change the dynamical evolution significantly provided that is sufficiently large and / is 
sufficiently small. All calculations are started with e = 0.9 and terminated when e = 0.2. 
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Figure 8. The dependencies of the orbital energy (the dashed curve), the rotational energy (the dotted curve) and the sum of 
the orbital and rotational energies (the solid curve) on time (7i„. All quantities are expressed in units of GMMp/ao- (Tin is set 
to 75. 

We start by discussing the case p = 2. In figure 7 we show the dependence of x = a/ao on 
time t. It is seen from this figure that the analytic solution (88) is in a very good agreement 
with the numerical solutions corresponding to (Tj„ = and 10. It also approximates rather 
closely the case of intermediate o-.j„ = 50. Results of calculation with high initial values of 
cr = 75 and 100 are in a good agreement with our analytic expression only for sufficiently 
large values of time t. In the case of a large (Tj„, initially, the semi-major axis (and eccen- 
tricity) grows with time II. Accordingly, the orbital energy also increases in the beginning of 
evolution. It turns out that the rotational energy of the primary is transfered to the orbital 
energy on a short time scale and this drives an increase of the semi-major axis. In figure 8 
we show the evolution of the orbital and rotational energies per unit of mass and their sum 
with time for the case = 75. The sum is gradually decreasing with time. That means 
that is it dissipated by the frictional processes. However, the orbital energy sharply increases 
and the rotational energy sharply decreases at time f < fi ~ 0.05. As it is seen from figure 
7 the moment of time t\ corresponds to the maximum of the curve x(i) with (7j„ = 75. In 
figure 9 we show the dependence of a on time for the cases with different (7j„. The solid curve 
corresponding to ai„ = shows a sharp growth of cr with time to the value a — Then the 
dynamical system is evolving in the cr = 8 resonance for a period of time t < ^2 ~ 0.23. After 

II Note that for all cases considered hereafter the orbital angular momentum is conserved in a good approximation and equation 
(73) is valid. 
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Figure 9. The same as Fig. 7 but the dependence of a on time t is shown. 

the time ^2 the system is evolving through a sequence of resonances with decreasing values 
of a. The solution corresponding to a"i„ = 10 is evolving in the a — 10 resonance from the 
very beginning. The late time evolution is similar to the previous case. The initial evolution 
of a for the case (jj^ = 50 looks like a monotonic decrease of a. In fact, as wc sec in the 
next figure 10 the system is evolving through a large number of resonances with decreasing 
resonance order n — 2a, but does not stay in any of these resonances for a long time. The 
curves corresponding to (7i„ = 75 and 100 are similar to the previous case with exception 
that these curves show that the initial value of a is maintained for a short initial period of 
time. This period of time corresponds to the stage of initial increase of the semi-major axis. 
In figure 10 we show the evolutionary tracks of our dynamical system on the plane {e,a). 
As it is seen from this figure the late time evolution proceeds through a sequence of the 
spin-orbit resonances with decreasing resonance order n. The dotted curve shows the ana- 
lytic dependence anie) corresponding to the case of the standard viscosity. It is important 
to point out that all curves corresponding to p = 2 have a > an for, practically, all values 
of eccentricity. Thus, a delayed response of convective motions to the tidal forcing can lead 
to a high state of rotation of the primary star. It is also very interesting to note that the 
dynamical system can evolve in a state of spin-orbit resonance even when the eccentricity is 
relatively small ~ 0.2. 

Finally, let us shortly discuss the case of p 7^ 2. In figures 11 and 12 we show the 
dependence of x on time i and the evolutionary tracks on the plane (e, cr) for p = 1.5 and 
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Figure 10. The same as Fig 7. and Fig. 8, but the dependence of cr(t) on eccentricity e{t) is shown. Note that the curves 
corresponding to (Jin = 50, 75 and 100 (the long-dashed, dot-dashed and dot-dot-dashed curves, respectively) arc close to each 
other and practically coincide at large values of e. The dotted curve represent the analytic solution cru{e) corresponding to the 
case of the standard Navier Stokes viscosity (see equation 75) 




Figure 11. The dependence of a; = a/ao on time t for different values of p. The solid, dashed and dot-dashed curves correspond 
to p = 2, 1.5, 2.5, respectively. (Tin is set to zero. The dotted curve represents the analytic solution (88). 



p = 2.5 in comparison with the previous case of p = 2. The initial value of the rotational 
parameter cTj„ is set to zero. From figure 11 it is seen that the case p = 1.5 {p = 2.5) evolves 
faster (slower) than the case p — 2 with respect to the time i toward the low eccentricity 
state. This effect is similar to that was found in the case of small p < 1 considered above 
(see figure 5). Figure 12 shows that the effect of evolution through a sequence of spin-orbit 
resonances remains valid for the different values of p > 1. The resonances are less prominent 
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Figure 12. The dependence of rT{t) on eccentricity e{t) is shown for p = 2 (the soUd curve), p = 1.5 (the dashed curve) 
and p = 2.5 (the dot-dashed curve). The dotted curve represent the analytic solution iTjy(e) corresponding to the case of the 
standard Navier Stokes viscosity (see equation 75). 

for the p — 1.5 case, and therefore the corresponding curve folfows the curve representing the 
standard Navier Stokes viscosity case on a large scale. On the other hand, the resonances 
corresponding to p = 2.5 case are very strong, and the corresponding dynamical system 
evolves in cr = 4 resonance until a small eccentricity e ~ 0.3 is reached. 

5 DISCUSSION AND CONCLUSIONS 

In this Paper we have developed a new self-consistent formalism for the calculation of the rate 
of energy and angular momentum transfer from an eccentric binary orbit to a fully convective 
primary planet or star. This is based on Fourier expansion in time of the tidal forcing 
potential and a normal mode expansion of the tidal response. To undertake calculations of 
orbital evolution, we have considered only the {I — 2) fundamental mode which dominates 
the response when the primary does not rotate, or when any potentially resonant inertial 
modes present when it rotates, have poor overlap with the forcing potential. 

We have assumed that the dissipation of the tidal response is due to an effective viscosity 
induced by convective turbulence. This is presumed to have an associated relaxation time 
or eddy turn over time which causes a delocalization of the dissipation process in time. It 
results in a weakening of the dissipation for tidal forcing at relative frequencies exceeding 
l/tc- We considered the case when the frequency dependence of the dissipative response is 
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oc 1/(1 + {uJm,ktcy), where ujm^k is the apparent frequency associated with the tidal forcing 
as viewed in the frame corotating with the primary in detaiL 

We use the fact that the orbital frequency is significantly smaller than the eigen- frequency 
of the {I — 2) f mode to introduce the equilibrium tide approximation. We note that, in 
general, this approximation is not equivalent to the classical assumption of a constant time 
lag between tidal forcing and response (eg. Alexander 1973, Hut 1981). However, we do show 
that they are equivalent only for the standard Navier Stokes viscosity law with instantaneous 
action of the viscosity (p = 0). 

We determine the orbital and rotational evolution of the primary numerically and an- 
alytically assuming that the angular momentum associated with the primary is much less 
than the orbital angular momentum for a range of values of p. We found that the evolution 
of the primary does not depend significantly on — D,tc, ^ being the orbital frequency, 
provided that ^ 1 and a scaled time t oc T^^t is used. 

However, the dependence on the parameter p was found to be significant. In the case 
of small p < 1 the evolution is similar to the case of the standard Navier Stokes viscosity 
{p — 0). When p > 1 and » 1, the system evolves through a sequence of specific spin- 
orbit corotation resonances with Q^/^ — ^/2, where Qj. is the rotation frequency and n is 
an integer. For p = 2we found an analytical expression for the evolution of semi-major axis 
given by equation (88). 

We stress that for p > 1 the primary rotation frequency at 'pseudo-synchronisation' is 
larger than obtained in the standard constant time lag approximation (p = 0) and that 
the primary can evolve in a state of spin-orbit resonance even when the eccentricity is small 
(e ~ 0.2). This may have observational consequences. In principal, one could use observations 
of the rotation of a sufficiently bright star in a close eccentric binary to test different models 
for response of convective turbulence to the tidal forcing. 

It should be possible to extend our formalism in several respects. For example, it can be 
generalised to take into account binaries with misaligned directions of orbital and rotational 
angular momenta (for the constant time lag case such a generalisation has already been 
made (eg. Alexander 1973, Hut 1981)). 

One can also extend the use of equations (34 - 35) to a situation where the contribution of 
other terms in the denominator (35) are non negligible compared with coq (as in fact occurs 
for sufficiently high harmonics). Taking into account a combination of several orthogonal 
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modes responding to the tidal forcing is also straightforward. This could generalise our 
formalism to include the lower order g modes in convectively stable stars. 

Instead of using the model frequency dependence of the viscous dissipation (equation 
(69)) one could hope to determine more realistic specifications from numerical simulation 
of convective media under a periodic forcing. Finally, one should use a realistic model of a 
planet or a star to find the response of the structure to the tidal heating and solve numerically 
the equations of evolution of the primary structure together with the equations of orbital 
evolution. 

Now let us discuss a possible restriction of our formalism due to the fact that the orbit 
may not be strictly periodic. Note that the binary period is changing with time as a result 
of energy transfer between the primary and the orbit. Let the change of the orbital period 
during one orbital period be APgrb and the corresponding change in the orbital frequency 
be An = -^^n. We see from equation (23) that the mode amplitude bm is expressed as 
a sum of harmonic terms, each of them being proportional to e*"^**, where uJk = kfl. When 
the number k is sufficiently large k > kmax, the change in the frequency Aujk — k\Sfl\ is 
larger than the difference between two neigboring frequencies equal to Jl. The expansion in 
Fourier series is clearly invalid for k > kmaxi where 



Now let us consider the resonance k — kr ^ cuq/D, which is found from the condition that the 
denominator R (equation (35)) is close to zero, and the corresponding bm,kr could be large. 
When kr > kmax, the corresponding term in scries (23) is not periodic, and the corresponding 
amplitude may be amplified with time. From the condition kr > kmax we find 



which, to within a numerical factor coincides with a condition for stochastic growth of the 
mode amplitude due to dynamic tides found by IP. When this condition is fulfilled and the 
energy transfer from the orbit to the mode is larger than the energy transfer associated with 
the equilibrium tide, our approach is not valid. 

Finally we comment that our results for the time to circularize an orbit starting from 
large eccentricity, tdrc, can be summarised for the various values of p in the compact form 



kr. 



"max 



APorb 



(90) 



(91) 



5MpU^ali 



(92) 



circ 




Using equation (71) for 7, this becomes 
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- ,^.J!l-ri!L.. (t.^/^Y . (93) 



37r(G'Mg)27o(l + q) 

The dissipation rate 70 is given in IP as 70 ~ 0.1/tc, inserting this into (93) and making 
some other reductions we obtain 

q'^al ( lOtct 



5 

^circ ^ r> 
671 



GM 



\l + q^ 



tjGM/al] . (94) 



The product of the quantities in square brackets in (94) can be estimated for a Jupiter mass 
protoplanet in the late stages of its evolution from data given in IP to be ~ 8. Similarly IP 
gives tc > lyr. Using these values and the representative value t — 0.2, we find for q — 10"^ 
that 

**'>-"(^)'(aI^)"('=V^)'- (-) 

This rather long time scale confirms the finding of IP on the basis of an impulsive 
treatment of orbits with very high eccentricity that equilibrium tides associated with the 
fundamental mode may not account for the circularization of the orbits of the recently 
discovered extra-solar planets with orbital periods of a few days (oq — 0.05AU) even when 
p = 0. 

Use of more plausible values, say, p > 1 would give circularization times well beyond a 
life time of the planetary systems ~ a few Gyrs. This may indicate that dynamic tides in 
the protoplanet (and possibly also the star for higher mass protoplanets ) and tidal inflation 
arising from dissipation in the protoplanet is important (see, eg. IP and references therein). 
Note, however, that the physics of convective turbulence is very poorly understood and the 
validity of the mixing length theory used is an issue. Also the contribution of the spectrum 
of inertial modes to the tide response may significantly shorten the circularization time scale 
(see eg. Ogilvie & Lin 2003). 
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APPENDIX A: APPROXIMATE ANALYTIC FORM FOR THE TIDAL 
FOURIER COEFFICIENTS 

As we discuss above the Fourier coefficients of the dimensionless tidal forcing amphtude (j)m,k 



Furthermore it is sufficient to consider only the case of m = 0, 2 and positive k, and for 
this case we approximately have am,k ~ Pm,k- Accordingly we find approximate expressions 
only for for a2,k, and q;o,a; hereafter. 

The integrals given in (Al) can be represented in a more useful form with respect to the 

eccentric anomaly ^ such that 



can be expressed in terms of two real quantities (j)jn,k — 2 



2\(^Tn,k Prn,k 



), where 




(Al) 




(A2) 



where 
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i?(e,cos(0) = 2e^ -1- 2ecos(0 + (2 - e'^) cos"^ {^) . (A3) 

The integrals (A2) cannot be evaluated analytically. Hence, in order to obtain an analytic 
representation, some reasonable approximations are needed. Often the integrands are devel- 
oped in powers of the eccentricity e and integrated term by term. However, the resulting 
series converge very slowly and are not convenient for analytic calculations involving orbits 
with even a moderate eccentricity e ~ 0.5. Here we use another approach based on the fact 
that the integrals are mainly determined by the region of integration near periastron ^ ~ 0. 
Accordingly, we expand the quantities R{e,cos{^)), D, (^ — esin(^)) in a power series in ^. 
Substituting these series in (A2) and truncating them, we get integrals which can be easily 
evaluated in principle. This approach gives good results in the limit of highly eccentric orbits 
e 1, and also for eccentricities as small as 0.2 and sufficiently large values of k. 

At first, let us discuss the evaluation of a2,k- In order to obtain an approximate value of 
Oi2,k we expand t(^) and D{^) in powers of ^ keeping terms up to third order. Then we may 
write 



/2(l-e)3^ x\ 

\ M ^ + y), (A4) 
D^{l-e){l + x'^), (A5) 
with X being a scaled form of ^ defined such that 



For an accurate approximation to a2,k, the quantity R{e, cos(^)) must be expanded up to at 
least fourth order so that 



i?(e, cos(O) - (1 - efil - ^(^x^ + ^ %^'\ '). (A7) 

Substituting equations (A4 - A7) into equation (A2) we get 

1 cos(A;r) , 2(2 + e) , 8 - e - 46^ , 

Finally, we need to express x in terms of r. For that, formally, one should solve the cubic 
equation ( A4) , but this would lead to significant technical difficulties. Instead we again make 
the assumption that the integral is mainly determined by the values of the integrand at small 
values of x and therefore r. This enables us to make use of the following approximate solution 
of equation (A4) which is correct to fifth order in these quantities and which also gives a 
linear relation between x and r for large x in the form 
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1 


2 


3 


4 







3.57-10-2 


2.05 • 10-2 


2.91 • 10-2 


1.34 • 10-3 


1.83 


• 10-5 


1 


6.76 ■ 10-2 


-1.68 ■ 10-2 


-1.13 • 10-2 


-9.74 • 10-4 


-2 • 


10-5 


2 


-9.97 • IQ-'^ 


-2.31 • 10-2 


2.115 ■ 10-5 


5.25 • 10-* 


2.18 


• 10-5 



Table Al. The coefficients a" determining (together with b") the values of integrals Ji given by equation (A12). Each row is 
labelled by the value of i, and each column by the value of n (see equation (A13)). 



where we adopt the re-scaled time 



Using equation (A9) to substitute for x in the integral (A8) we get 

2 / 2(2 + e)^ 8-e-4e2 \ 
«2, - ^ Uo - + ^^J^) > (All) 



where 



~2\5-2i~2i/i I 2~2\2 



Ji 



POO M j_ -7- r i + -r 

/ df^^^^ ^ -t ' ' cos{k\f). (A12) 



Here, again making use of the assumption that the integral is mainly determined in a small 
interval near f = 0, and therefore to a good approximation independent of the interval it is 
taken over, as long as that is big enough, we have changed the limits of integration to ±oo. 

Note that we have also neglected terms proportional to in the denumerator of the 
integrals (A12). These terms may easily be taken into account in principle but do not change 
significantly the results. The integrals (A12)) can be easily evaluated using the theory of 
functions of a complex variable. They are determined by two poles of fifth order in the 



complex r upper half-plane located at r = (3 =F \/(6)). The Jj so found can be represented 
in the form 

?i=4 n=4 

J. = 6-2/+^^^ ali^kf + e~y-^^ Y ^r(A^)", (A13) 

n=0 n=0 

where y± — ^(3 =F \f{^)), and the coefficients af and 6" are given in Tables Al and A2 
respectively. 
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1 


2 


3 


4 





-4.77 • 10-2 


-5.17-10-2 


-6.12 • 10-3 


1.33- 10-2 


5.64 • 10-3 


1 


-9.25 ■ 10-=* 


4.54 ■ 10-2 


7.04 ■ 10-2 


2.29 ■ 10-2 


-1.08 ■ 10-2 


2 


5.84 ■ 10-2 


1.52 • 10-2 


6.405 ■ 10-2 


-1.36 • 10-1 


2.05 ■ 10-2 



Table A2. Same as table Al but for the coefficients 61* 



n 


1 


2 n 


1 


2 


c„ 0.499 


0.141 


6.38 ■ 10-3 dn 


-3.81 • 10-2 0.141 


0.199 



Table A3. Coefficients Cn and dn determining the value of the integral (A15). 

We evaluate ao^k by exactly the same method. Doing this we get 
where 

7=rdfii±l2^^. (A15) 

The integral (A15) is determined by the behaviour at two poles of third order located at 
f = iy± in the upper complex half-plane and can be expressed as 

n=2 n=2 

/ = e-^+^*^ ^ni^^T + e""-'' E (A16) 

n=0 n=0 

where the coefficients and dn are given in table A3 . 

This paper has been produced using the Royal Astronomical Society/Blackwell Science l^T^ 
style file. 
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